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Abstract. We introduce a new aproximation scheme for the periodic Anderson model 
(PAM) . The modified alloy approximation represents an optimum alloy approximation 
for the strong coupling limit, which can be solved within the CPA-formalism. Zero- 
temperature and finite-temperature phase diagrams are presented for the PAM in 
the intermediate-valence regime. The diversity of magnetic properties accessible by 
variation of the system parameters can be studied by means of quasiparticle densities 
of states: The conduction band couples either ferro- or antiferromagneticaly to the /- 
levels. A finite hybridization is a necessary precondition for ferromagnetism. However, 
too strong hybridization generally suppresses ferromagnetism, but can for certain 
system parameters also lead to a semi-metallic state with unusual magnetic properties. 
By comparing with the spectral density approximation, the influence of quasiparticle 
damping can be examined. 



PACS numbers: 75.30.Mb 71.28.-fd 71.10.Fd 75.20.Hr 
1. Introduction 

Correlated electron systems have come to occupy the centre stage in both theory and 
experiment in condensed matter physics. There is a general consensus that strong 
electron correlations play a decisive role in a variety of phenomena such as magnetism, 
heavy fermions, high-temperature superconductivity, colossal magneto-resistance etc. 
The real systems that display these phenomena are transition metals (3d), alloys 
and compounds as well as rare earth (4f) systems (metals, insulators), in particular 
cuprates, manganites, Ce compounds etc.. There is intense activity both theoretically 
and experimentally to isolate the essential interactions responsible for these phenomena. 
In this process, several model systems are being studied which are expected to at least 
qualitatively mimic the actual systems in some region or the other of the parameter 
space. One of the most widely used models which brings out the role of electron 
correlations while also taking into account the interplay of two different types of 
electrons, one highly localized and the other quasi free, is the well-known Anderson 
model fl]]. One distinguishes the "single impurity Anderson model" (SIAM), if the 
system of uncorrelated conduction electrons hybridizes with a single localized state, and 
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the "periodic Anderson model" (PAM), if the hybridization takes place with a periodic 
lattice of localized states. As far as the SIAM is concerned many exact results have 
been obtained, e. g. by using Bethe ansatz or renormalization group theory . Very 
recently the SIAM has gained a new upsurge in connection with the "dynamical mean 
field theory" 0, ^, ^ that exploits the fact that in infinite lattice dimensions theoretical 
lattice models like the Hubbard model can be mapped onto single-impurity models such 
as the SIAM. Numerically essentially exact methods such as the exact diagonalization of 
small systems, quantum Monte Carlo calculations or numerical renormalization group 
theory have provided further insight into the physics of the model. However, since 
these numerical methods are bound to certain limitations, reliable analytical approaches 
remain to be required p. 

The situation is less satisfactory in the case of the PAM. Even though there exists a 
number of approximate schemes for different parameter constellations, referring to either 
of the so-called heavy fermion, intermediate valence and Kondo regimes, there is still 
need for further improvement and extension. The matter acquires additional urgency 
in view of the fact that for the parameter space corresponding to strongly correlated 
low-lying localized states the PAM can be mapped |^ on to the Kondo lattice model 
||T0| which is of great current interest, e. g. with respect to the extraordinary physical 
properties of the manganites . 

Recently we proposed a new approximation scheme which is based on a mapping 
of the PAM on to an effective Hubbard model. The parameters of the Hubbard model 
are such that they correspond to the strong coupling limit. Therefore we could exploit 
a reliable approximate theory, trustworthy in the strong coupling limit. Such a theory 
is the "spectral density approach" (SDA) [|1^, [1^. Its main advantages are the 
physically simple concept and the non-perturbative character. In this scheme we have 
studied the magnetic T = 0-phase diagram of the PAM as well as its finite temperature 
magnetic properties. The SDA, however, shows up a principal limitation concerning 
quasiparticle damping. By ansatz the SDA self-energy is a real quantity, thus neglecting 
from the very beginning effects due to finite lifetimes of the quasiparticles. How 
quasiparticle damping influences the magnetic stability in the PAM is an interesting 
and important question having been left open by our previous theory |]12[. It is the 
main impetus for the present study to close this gap. 

A concept ionally simple method that is able to provide complex self-energes is the 
"coherent potential approximation" (CPA) [|1^. In order to perform CPA on the PAM, 
we have to fix up an alloy analogy. This requires the determination of the energy levels 
and the concentrations of the components of the fictitious alloy. Within the conventional 
alloy analogy treatment (AA) due to Hubbard [|1^] these are taken by referring to the 
atomic limit. For an application to the PAM, see e. g. |]TB|, r^|. However, this choice is by 



no means predetermined. On the contrary it can be shown that it is indeed not the best 
ansatz. One knows that CPA becomes an exact procedure in infinite lattice dimensions 
{d = oo), where the inherent single-site aspect of CPA is rigorous However, the 
CPA-solution of the AA for the d = oo-Hubbard model violates the exactly known 
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strong coupling behaviour |^ and also does not reproduce the weak coupling results of 



second order perturbation theory p2| , |23[] . Furthermore it contradicts exact high-energy 
expansions |2|]. So we have to conclude that the conventional AA, which starts from the 
atomic limit solution, is not the most convenient alloy analogy. In this paper we therefore 
derive the "modified alloy analogy" (MAA), which was recently introduced for the 



Hubbard model P3|, for the PAM. This method substantially improves the AA ansatz 
by deriving the proper alloy analogy from exact high energy expansions of the single- 
electron Green function and the self-energy, respectively. This procedure guarantees 
the correct strong coupling behaviour relevant for the PAM. On the way of deriving 
the MAA, the above-mentioned SDA can easily be re-derived without the necessity of 



introducing the effective Hubbard model as was done beforefT^. It turns out that the 
"atomic levels" of fictitious alloy constituents are nothing else but the SDA-quasiparticle 
energies in the zero-bandwidth limit. Similarly the "concentrations" agree with the 
SDA-spectral weights in this limit. The MAA therewith promises to keep the advantages 
of the SDA while simultaneously improving the method by incorporating quasiparticle 
damping. The SDA-energies and spectral weights contain non-trivial thermodynamic 
expectation values which have to be determined self-consistently. By them the itineracy 
of the — (T-electrons, which define the fictitious alloy for the propagating cx-electrons, 
comes into play, at least to a certain degree. The neglect of this itineracy is a well- 
known shortcoming of the conventional AA. 

The paper is organized as follows. In the next section the PAM and its many- 
body problem is formulated. An approximate solution by use of the MAA scheme is 
then developed in section 3. Section 4 is devoted to a presentation and discussion of 
the results concerning the magnetic properties of the PAM. The paper ends with some 
concluding remarks. 

2. Model Hamiltonian and Its Many-Body Problem 

Starting point is the periodic Anderson Hamiltonian 

H = - fi)sls,, + ^(ef - f^)fU, + (1) 

ija i,cr 

+ V Y.^fls.. + .1/..) + \UY^ n!£r^^ (2) 

Sio- (/to-)and {fla) ^^6, respectively, the annihilation and the creation operators for 
an electron in a non-degenerate conduction band (localized /-state), and n\l) = fj^fia 
is the spin-dependent occupation number operator for the f-state. The index i refers to 
the respective lattice site, a = {t, i} is the spin projection. The hopping integral T^, 

k 

describes the propagation of a band electron from site Rj to site e(fc) is the free 
Bloch energy, while Cf denotes the position of the non-degenerate /-level. We choose 
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the energy zero to coincide with the centre of gravity of the unperturbed conduction 
band: 

T.. = ^$^6(fc) = (4) 
fc 

U is the intra-atomic Coulomb repulsion between /-electrons. The hybridization V is 
taken as a real and fc- independent local matrix element, n is the chemical potential. 
Let us start with the retarded single-s(/) electron Zubarev Green's function 

GiP^E) = {{U; fl))- Gi;l{E) = si)) (5) 
fc 

The equations of motion are easily derived and formally solved, 

"""^ ^ (E - (Cf - /i) - ^UmE - (e(fc) -f^))-V^ ^ ^ 

GkliE) = y2 ^ .px (8) 

E-{ei- fi}- E-{e{k)-i,) ~ ^ka[E) 

where the self-energy has been introduced via 

^kAE)G^£iE) = f/1 $^((4/,-./p+fc-,.; flj (9) 

Obviously the explicitly problem is solved as soon as we have found a solution for the 
self-energy. The /- and s-quasiparticle densities of states (QDOS) are given by 

Pi'\E) = Yl ^^i^ - + iO"") (10) 

fc 

Pi'^ (E) = ^^i'i - + iO^) (11) 

fc 

(s f) 

and determine the spin-dependent average occupation number Ua , which we need to 
construct the magnetic phase diagram of the PAM: 

/ + CO 
dEf4E)pi'\E) (12) 
-oo 
r+oo 

= iflU) = / dEf4E)pif\E) (13) 



3. Modified alloy analogy 

A standard method to solve many-body problems as that posed by the PAM 
Hamiltonian (|l|) is the CPA |T^. As a single-site approximation the resulting /-electron 



self-energy will be local, i. e. fc-independent. If we apply this method here we have to 
solve the following equation self-consistently 

Epa — T.o-{E) — Cf 

t^i '"l-lG."£(B)(S-E,(B)-e,) 
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The solution of (0) needs to fix the alloy analogy, i. e. the "atomic levels" Ep^ and the 
"concentrations" Xp^ of the n constituents of the fictitious alloy. Since the dominant 
features of the correlated /-level in the PAM are the two charge excitations separated 
by U, the number of "alloyed" components n is set to two. The conventional alloy 
analogy (A A) uses the = limit ("atomic" limit) of the PAM to determine Ep„ and 

Eir' = eu = (16) 

^ = e, + U- x(f ) = 1 - nL^] 

The result violates the weak as well as the strong coupling behaviour and as for the 
Hubbard model prohibits spontaneous magnetism [ll8| , |26| , |2^ , In particular the 
strong coupling behaviour appears to be crucial with respect to ferromagnetism. On 
the other hand, the choice (|T^) is not at all predetermined. We propose another way to 
find out the "best alloy analogy". Correct strong coupling behaviour is guaranteed by 
fulfilling the high-energy expansion of relevant Green functions and self-energies 



Decisive ingredients for proper high-energy expansions are the local spectral moments 

M^J'^ = ]^Y.Mu^-^ ^ = 0,1,2,... 
fc 

M^-j = i 1 dEE-Su^E) (17) 
Sko-{E) is the /-electron spectral density: 

SUE) = --QG^£{E + iO+) (18) 

TT 

The moments can be calculated independently of the required spectral density: 

MjrJ = {[ [... [h.,H].,...,HU fLU) (19) 

n-fold commutator 

[...,...]_ denotes the commutator and [...,...]+ the anti-commutator. For the Green 
function (|^) we can write 

fdE'^ = n±M^ (20) 

For solving the CPA-equation (|Tip we only need the local spectral moments which are 
found to be 

m(°) = 1 

M^i) =€{ + Un^H 

M^2) = ^ 2ef[/nL^] + U^n^H + 1/^ (21) 
) = ef + ^ejUri^H + [/^ef(2nL^] + n^lf) + U'^n^H + 
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From d^) we get the corresponding expansion for the self-energy: 

oo ^(n) 

S'^'^(^) = E#' ^'^(E) = j;^J2^>.AE) (22) 

n=0 k 

with the local coefficients: 

cm = Un^H 

C^^) = U^n^ll{l-ri^ll) (23) 
Cf) = f/2^(/)(i _ rSll)[B_^ + f/(l - n^H)) 

Surprisingly the hybridization V does not explicitely appear in the Ca^\ The 
contributions via the moments (|21| ) are exactly cancelled by those from the term 

By use of (^) and p^ ) we now expand the CPA-equation ( p!^ with respect to 
powers of Comparison of the coefficients of the terms up to n = 3 yields the 
following set of four equations for the four unknown quantities -E(i,2)o- and X{i^2)a'- 

2 

^ ^ Xpcr 1 
p=l 
2 

Y,Xp^{Ep. - e{) =Un^ll (24) 
p=i 

2 



p=i 

2 

^Xp.(Ep. - e,f = U'n^Il + [/2nL^](l - nL^])(fi_. - Cf) 



p=i 

We now use eqs. (^) to fix the "optimum" alloy analogy. After simple manipulations 

we get 

Ei,2. = \ [B-a + f/ + ef ± \J (5_, + U- e,f + AUn^ll{e, - B_„)] (25) 
_ - ef - t/nL^] _ ^ 

a^lcr — ^ — J- — 3^20- 

The decisive point for the ongoing procedure is the "higher" correlation function -B-^- 
n^Hil - - ef) = l^(/E^.-.(2aS^ - 1)) (26) 

It has been elaborated in previous works [|1^, |1^, |4| that its analogue has a decisive 
influence on the stability at spontaneous magnetism in the Hubbard model. We believe, 
it plays a similarly important role for ferromagnetism in the PAM|jT2|, In spite 
of the fact that it is a "higher" correlation function it can rigorously be expressed by 



single-electron terms [0 



n^llil - n^ll){B_^ - ef) = -—^ J dEf_{E) [jj^AE) - ij (27) 



Ferromagnetism in the Periodic Anderson Model 

x(^{E-{e,-fi)-J:^{E))G\2iE)-h 

With the fictitious alloy (pSD, we enter the CPA-equation ([I^). The theory is now 
complete. The equations (|l|), (|15D, (D, (^) and (|13D, together with (H) build a closed 
system of equations which can be solved self-consistently for the self-energy ^^-(i?). 
Note that the only fc-dependence comes into play by (|^) via the formal solution of the 
equation of motion. The fc-dependence is therefore strictly an e(fc)-dependence, so that 
all k-summations can be replaced by simpler energy-integration over the "free" Bloch 
density of states po{E) = ^/Nj^k^i^ ~ ^(^)) which has to be considered as a model 
parameter: 

G^l (E) = h dx (28) 

J-oo i?-(ef-M)-^4^-S.(E) 
Let us comment on the modified alloy analogy (pSf) and its differences to the 



conventional AA (|T6D: Within the AA, the a-f electron is propagating through a 
fictitious alloy where one component is represented by lattice sites with no — a-electron 
present, and the other component by lattices sites which are occupied by one —a 
electron. This approach therefore completely neglects the hybridization between the 
/-levels and the conduction band. For calculating the self-energy, the —a /-electrons 
are frozen, any exchange with the conduction band is eliminated. This excludes the 
possiblity of magnetic order in the system. How can this drawback be circumvented? 
It is clear that the PAM will always be considered to belong to the strong coupling 
regime (large U). Two well separated excitation peaks will therefore be expected in 
the /-electron quasiparticle density of states. The above- described ansatz of fitting the 
positions and weights of these to an exact -^-expansion of the Green function seems 
the most plausible one since these charge excitation peaks themselves are high-energy 
features. It is noteworthy that all thereby induced differences between the MAA and the 
AA are due to the hybridization between /- and s-levels. It can be seen in equation (^) 
that in the limit of vanishing hybridization (V^ = 0), B^^^ = Cf holds and equations (^) 
reduce to the conventional alloy analogy (|16]). However, for finite hybridization, B_„ 
and n^Il are to be determined self-consistently by (pTj) and (|13|), respectively, possibly 
providing the alloy energies and concentrations with an explicit spin-dependence. In 
each step of the CPA-iteration and may change and so does the underlying 
fictitious alloy. Equation (^) makes clear that the inclusion of accounts to a certain 
degree for the hybridization of —a- /-electrons with the conduction band and therefore 
also to their effective itineracy which is completely neglected in the conventional alloy 
analogy (|16]). 

Before further discussing the MAA, let us remind the reader of the already 
mentioned spectral density approximation (SDA). In a previous work, this theory was 
derived via mapping of the PAM onto an effective medium Hubbard model |T^ . Within 
this effective Hubbard model, the SDA is justified in the strong coupling limit where the 
obtained positions and spectral weights of the lower and upper Hubbard band coincide 
with exact results up to the order ^0]. At this point it is noteworthy that exactly the 
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same results as in reference can also be obtained by adapting the concept of the SDA 
directly onto the PAM which is possible and straightforward from the above-presented 
results. Motivated by the solution of the "atomic" limit {V = 0) we make the following 
ansatz for the self-energy: 

= (29) 

-c/ — 

The coefficients of this ansatz can now be fitted in such a way that the high energy 
expansion of the self-energy, (^) with the coefficients (23) is fulfilled. One readily 
arrives at 

^(SDA) ^ Un^IliE - B^^ - e,) 

E-B_^-et-U{l-n^Il) 
which together with equations (0), ( pTj) and (|T3|) solves the problem. It is 
straightforward to show that this is identical to the approximation proposed in reference 
||12|| . Therefore, the discussion of the advantages and disadvantages of the method found 
therein remains valid. The main disadvantage was clearly the neglection of quasiparticle 
damping in the ansatz (p9|) . 

Now turning back to the MAA, its benefits are immediately clear. While 
reproducing the high-energy expansion up to the same order as the SDA, it additionally 
incorporates quasiparticle damping via the CPA formalism. And contrary to the 
conventional alloy analogy, the freedom of defining the " alloyed components" is used 
to ensure the correct high-energy expansion (|22|) . 

We conclude that the essentials of the qualitatively convincing SDA used in 
our previous paper [1^ are incorporated in the MAA and completed by a proper 
quasiparticle damping. So the MAA represents a systematic improvement and extension 
of the SDA. Comparison of results for MAA and SDA will allow to inspect very directly 
the influence of quasiparticle damping on magnetic stability in the PAM. 



4. Results and Discussion 



The presented theory for the periodic Anderson model has been evaluated for a system 
characterized by the following parameters: The free conduction band density of states 
is chosen to be semi-elliptic with a width of W = leV centred a.t E = 0. The /-level 
is determined by its distance (cf) to the conduction band centre and the intraatomic 
Coulomb interaction U. The latter is 4eV throughout the paper. That means we 
consider the PAM in the strong coupling regime. Furthermore, we are mainly interested 
in the intermediate valence regime, i. e. Cf is presumed to be located within the Bloch- 
band region {—0.5eV < Cf < +0.5eV). It turns out that the physics of the PAM strongly 
depends on the total particle density 

^(tot) _ ^(^W +4^)); < < 4 (31) 

We therefore present results for different n^*"*-*, but so that the upper charge excitation 
at Cf + f/ ("upper Hubbard band") remains in any case unoccupied (n*^*"*) < 3). Another 
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decisive variable is of course the temperature T given in Kelvin. The evaluation of our 
theory has been done for a translational symmetric lattice, antiferromagnetic ordering 
is not considered. We rather concentrate ourselves on the possibility and the conditions 
of spontaneous ferromagnetism. 




Figure 1. Magnetic phase diagram in the n^*"*^ - ef plane for V — 0.1, U — 4 
and T — O.f— and s— magnetizations are parallel to the left of the dotted line and 
anti-parallel to the right. 



The magnetic phase diagram in terms of n*^*°*^ and Cf for a small hybridization 
V = O.leV is plotted in figure 1. For a given Cf within the Bloch band region 
ferromagnetism becomes possible for particle densities in between a lower and an 
upper critical value. These critical values shift to higher numbers with increasing 
/-level position. Ferromagnetism is of course basically provoked by the correlated /- 
electrons. Except for the hybridization with the conduction electrons they are described 
by the zero-bandwidth Hubbard model. The latter special case, however, forbids 
ferromagnetism. So a finite hybridization V is needed. On the other hand, one knows 
that ferromagnetic spin order in the Hubbard model is bound to further conditions: 
The band occupation and the Coulomb correlation U/W must exceed critical values 



which are different for different lattice structures [|T^ [T5|, 0. Furthermore, the 
free band should have an assymetric density of states|2T|. SDA and MAA both do not 
allow ferromagnetism in the "pure" Hubbard model for highly symmetric free densities 
of states such as the semi-elliptic one irrespective of the correlation strength U/W. In 
the PAM with finite hybridization V, however, all these conditions can be met: The 
hybridization creates a finite width of the /-dispersion, and the resulting effective /-band 
turns out to be strongly asymetric, therewith allowing for a ferromagnetic ground state. 
For a given ef the lower critical density in the phase diagram (figure 1) is the analogue 
to the critical particle density in the Hubbard model [0. At the upper boundary the 



lower /-charge excitation is more or less filled corresponding to a half-filled band in the 
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Band filling n Band filling n Band filling n' 



Figure 2. /— (solid line) and s— (dotted line) magnetizations as a function of n^*"*-* 
for various values of Cf. The s— magnetization is multiplied by a factor of 5 for better 
visibility. Other parameters are as in figure f . 



strongly coupled Hubbard model for which antiferromagnetism is to be expected]^, [T^. 

The upper boundary in the phase diagram for a given Cf corresponds always to a 
first order transition, the lower one to a second order transition. This can be seen in 
figure 2 where the spontaneous T = 0-magnetization is plotted as function of the total 
particle density n^^°^\ The position of the /-level leads to strikingly different behaviour 
of the magnetization. 

With Cf in the lower half of the conduction band (left column in figure 2) 
the /-magnetization increases from the lower boundary very rapidly into saturation 
in order to perform the mentioned discontinuous transition into the paramagnetic 
phase at the upper boundary. Interesting is the induced s-polarization starting 
positive (ferromagnetic s-f coupling) for weak electron densities and changing its sign 
(antiferromagnetic s-f coupling) at about = 1.2. According to the Schrieffer- 

Wolff transformation at first glance an antiferromagnetic coupling between s-and /- 
electrons would be expected. However, the transformation is bound to the Kondo-limit 
(7^{/) = X) and is not necessarily conclusive for the intermediate valence region which 
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-0.5 0.5 -0.5 0.5 

Energy E Energy E 

Figure 3. /— and s— quasiparticle densities of states for various values of n*^*°'^ and 
Cf = —0.3. Full lines for spin-up and dotted lines for spin-down. Other parameters are 
as in figure 1. Thin vertical lines show the position of chemical potential. 



is investigated here. The quasiparticle density of states (Q-DOS), plotted in figure 3 
for Cf = —O.SeV and various electron densities, gives some evidence how to understand 
the sign change of the s-electron polarization. It is instructive to decompose the Q- 
DOS into s- and /-parts, although this decomposition is artificial because of the finite 
hybridization. The latter takes care for the fact that s- and /-partial-Q-DOS occupy 
exactly the same energy regions but of course with different weights. In the j-parts of 
both spectra the hybridization gap around Cf = —O.SeV is clearly visible. For practically 
all densities n^^°^\ exhibited in figure 3, the /-system is saturated, i. e. there are no 



down-spin / electrons. It is known from the Hubbard model |T^, |I5| and will be shown 
in figure 13 for the PAM, that then the damping of up-spin quasiparticles is in such a 
case negligible compared to that of down-spin quasiparticles. Fine structures like the 
hybridization gap are to be seen in the up-spectrum but not in the down-spectrum. 



The comparison with the SDA results in ref. [|I2[ indeed confirms the interpretation 
that quasiparticle damping closes the hybridization gap in the |-spectrum. For the 
SDA neglects such damping effects, consequently the SDA-|-spectrum, too, exhibits a 
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Energy E Energy E 

Figure 4. Same as figure 3 but ef = —0.1. 

gap. There is a slight exchange shift in the s-part due to hybridization. As long as 
the chemical potential /i is below the hybridization gap the system contains more up- 
than down-spin s-electrons, the s-f coupling is ferromagnetic. When fi shifts above the 
gap (in figure 3 in between n^^°^^ = 1.10 and n*^*"*-* = 1.25), which exists only in the 

(s) (s) 

up-spectrum, then the increase in distinctly exceeds that of , the s-f coupling 
becomes antiferromagnetic. That explains the magnetization behaviour in figure 2 (left 
column) as a density of states effect. Note, however, that the absolute value of the 
induced s-polarization is always rather small. 

A qualitatively different situation is observed for the magnetization when Cf is near 
to the centre of the s-band (middle column in figure 2). A reentrant course appears, 
obviously not connected with the sign change of the s-polarization. It manifests itself 
in the phase diagram of figure 1 by the "oscillating" phase boundaries. Again a look at 
the respective Q-DOS is quite instructive, figure 4 shows the example Cf = —O.leV. For 
such /-level position the hybridization gap splits the /-dominated part of the spectrum 
rather symmetrically into two almost equally weighted peaks. That holds in particular 
in the f-spectrum of the ferromagnetic phase (n^*"*^ = 1.35, 1.45, 1.75, 1.85 in figure 4). 
The chemical potential /i is located in one of the peaks, a situation which according to the 
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Figure 5. Same as figure 3 but ef = 0.3. 
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simple Stoner criterion favours the appearance of ferromagnetism. As already discussed 
for the case of figure 3 (left column in figure 2) in the ordered phase f-quasiparticles are 
rather long-living while the |-particles are strongly damped. Again we conclude that this 
is the reason why the hybridization gap exists only in the |-spectrum. The spontaneous 
ferromagnetism disappears when /i enters the hybridization gap in accordance with the 
Stoner criterion. The quasiparticles are stronger damped in the paramagnetic phase 
than the f-particles in the ordered phase, so that the hybridization gap is more or less 
covered by the "smeared-out" quasiparticle peaks. 

When the bare /-level Cf has been shifted into the upper half of the conduction band 
the polarization type changes once more. The coupling between s- and /-electrons is 
always antiferromagnetic (right column in figure 2). The Q-DOS in figure 5, concerning 
the case Cf = 0.3eV, reveals that ferromagnetic order sets in only when fi is above the 
hybridization gap. The consequence is that there are more s-electrons with spin | than 
spin |. 

The phase diagram in figure 1 is calculated for a fixed hybridization V — O.leV. 
It is clear that V has a decisive influence on the extension of the ferromagnetic phase 
in the ef-n^*°*^ plane. Figure 6 demonstrates that the higher the hybridization V the 
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Band filling n Band filling n 

Figure 6. /- and s- magnetizations as a function of for different values of 

hybridization strength V with Cf = —0.3. Other parameters are as in figure 1. 



smaller the ferromagnetic region. On the other hand, the induced polarization of the 
s-electrons increases with growing hybridization. Antiparallel coupling appears only for 
weak V. Strong electron fluctuations between s-band and /-level generally diminish 
magnetic stability but enhance the polarization of the s-band. However, the same 
hybridization is exclusively responsible for the ferromagnetic order. V = does not 
permit ferromagnetism. 




0.5 1 1.5 2 2.5 



Band filling n 

Figure 7. Magnetic phase diagram for MAA and SDA with V = 0.2, [/ = 4 and T = 0. 
/— and s— magnetizations are parallel to the left of the dotted line and anti-parallel 
to the right in the SDA. In the MAA /— and s— magnetizations are ant-parallel only 
when > 2. 



It is an interesting point to find out to what an extent quasiparticle damping may 
influence the possibility and the stability of ferromagnetic order. This can best be done 
by a comparison of the results found by SDA[Q and by MAA. Both methods are based 
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on the same physical ideas, MAA can be classified as " SDA plus quasiparticle damping" . 
Figure 7 shows the magnetic phase diagram, derived within, respectively, SDA and MAA 
for the same set of model parameters. The Cf-n^*"*-* area of the ferromagnetic phase is 
distinctly restricted in the MAA because of the quasiparticle damping, contrary to the 
SDA-result which is free of damping effects. For both methods we find that for lower 
particle densities a ferromagnetic, for higher densities an antiferromagnetic coupling of 
the s-electron spins to the local /-moments takes place. 
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Figure 8. /— (solid line) and s~ (dotted line) magnetizations as a function of 
hybridization strength V for different n'*"*^ values for Cf = —0.3 at T = 0. The 
inset shows the s~ magnetization for = 1.2 in the region of small V. 



Up to now we discussed mainly the Cf- and the n(*°*^-dependence of the physical 
properties of the periodic Anderson model. Let us now examine the influence of V in 
detail. Some results are plotted in flgures 8, 9 and 10 for Cf = —0.3eV and several 
particle densities n^^°^\ The behaviour is not at all unique. For rather 
weak hybridization V is sufficient to destroy the ferromagnetism (figure 8). The Curie 
temperature Tc runs through a maximum (figure 9) pointing therewith to the fact that V 
provokes two competing effects. It broadens the /-level, thus creating the precondition 
for a magnetic order (no ferromagnetism in the zero-bandwidth Hubbard model [pO|]!). 
On the other hand, increasing s-f fluctuations must damage the ferromagnetic order 
because the empty /-level does not carry a magnetic moment. It is remarkable that Tc 
goes smoothly to zero (flgure 9) although the T = 0-magnetization performs a flrst-order 
transition (flgure 8). 

A completely different ^-dependence shows up for the band occupations n*^*"*) = 1.0 
and 1.2. In these cases large hybridization strength does not destroy the ferromagnetism, 
but rather enhances Tc. By closely examining flgure 8, one recognizes a "kink" or a 
minimum in the respective magnetization curves, for = 1.0 at 1/ ^ 0.45 and 

for n(*°*) = 1.2 at V ^ 0.09. By comparing with the densities of states (flgure 10), 
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Figure 9. Curie Temperature as a function of V for different values of n*^*°'^ with, 
ef = -0.3. 



we note that this value approximately coincides with the l^-value where \i enters the 
hybridization gap in the |-Q-DOS. I. e. the system becomes a semimetal where only 
the |-electrons contribute to the electrical current. One might speculate that this has 
something to do with the increasing T^. However, we cannot rule out the possiblity that 
this is an effect of our approximation rather than an inherent property of the PAM. 

Let us finally discuss the temperature-dependence of the magnetic properties as, 
e. g., the spontaneous magnetization (figure 11). First order as well as second order 
transitions appear. Whether first order transitions are artefacts of our approximate 
procedure or intrinsic properties of the PAM is not clear. A similar situation is found 
when applying the MAA to the Hubbard model [|^. In case of a second-order transition 
the /-magnetization behaves like a Brillouin function, while the corresponding induced 
s-polarization very often shows remarkable deviations. The examples plotted in figure 
11 demonstrate again that the s-j coupling may be ferromagnetic or antiferromagnet 
depending on the band filling. 

A key-quantity of ferromagnetism is the Curie-temperature Tc, which is of course 
decisively influenced by the bandfilling (figure 12). The 72^*°*-' -dependence is very much 
more regular than the V^-dependence exhibited in figure 9. The reentrant behaviour for 
Cf = —^.\eV corresponds to that of the T = 0-moment in figure 2 (middle column). 
The calculated Tc-values are of a realistic order of magnitude. The transition into the 
paramagnetic phase {T^ = 0) seems to be always continuous even if the break-down 
of the T = 0-moment is discontinuous (figure 2). The temperature-dependence of the 
magnetization is due to a respective behaviour of the Q-DOS (figure 13). There is a 
distinct spin asymmetry in the lower /-like peak (" lower Hubbard band" ) which causes 
the spontaneous magnetic moment. The up-spin part is characterized by a hybridization 
gap, which is not visible in the down-spin spectrum (cf. figures 3, 4 and 5). As already 
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Figure 10. /- and s- quasiparticle densities of states for various values of V with 
Cf = —0.3 and = 1.2 at T = 0. Full lines for spin-up and dotted lines for 

spin-down. Thin vertical lines show the position of chemical potential. 

mentioned, this is due to quasiparticle damping, which for low temperatures is very much 
stronger for |- than for f-quasiparticles (cf. figure 13). With increasing temperature 
the damping of the t-pai'ticles, is getting larger and that of the J, particles smaller. 
For increasing temperatures, a dip develops in the spin-| Q-DOS, which eventually 
merges with the spin-| hybridization gap at T = T^. The spin asymmetry is then 
removed. The induced spin polarization of the conduction band is always very weak 
so that the assumption that the collective order is based on an RKKY-like coupling 
via the polarized conduction electrons appears unlikely. It is rather to believe that the 
observed ferromagnetism is due to strong electron correlations in the narrow "/-band" 
as it happens in the single-band Hubbard model [^, The coupling to the s-band 
via V takes care for the finite width of the original /-level being in this sense the basic 
precondition for ferromagnetic order in the periodic Anderson model. 
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temperature for various Cf and as mentioned in the figure, s— mgnetization 
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Figure 12. Curie temperature Tc as a function of band filling n^*"*-' for different values 
of Cf. The hybridization strength V = 0.1 for the left and V — 0.2 for the right of 
figure. 



5. Conclusion 

A "modified alloy analogy" (MAA) p5|, previously introduced as an approach to 



the single-band Hubbard model 0] has been applied to the periodic Anderson model 
(PAM). A high-energy expansion of relevant Green functions as well as a corresponding 
expansion of the determination equation for the CPA-self-energy have been used to 
find the optimum alloy analogy for the PAM. This alloy analogy is then used to solve 
the PAM-many body problem within the CPA. By construction the MAA represents 
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Figure 13. /- and s- quasiparticle densities of states at various temperatures. ef=- 
0.3, F=0.1 and =1.1. Full lines for spin-up and dotted lines for spin-down. Thin 

vertical lines show the position of chemical potential. Additionally, the imaginary part 
of the self-enery is plotted as thin lines in the right column. For better visibility, it is 
multiplied by —100. 



a strong-coupling approach, therewith most probably suitable to describe spontaneous 
ferromagnetism in the PAM. It can be considered as an extension and improvement of the 

|], mainly by inclusion of quasiparticle 



'spectral density approach" (SDA) [|13], [Tj, |T5 



damping. It cannot reproduce, however, the low-energy features of the PAM (Kondo 
resonance) im, being certainly not so decisive with respect to magnetic stability. It 
incorporates, however, important higher correlation functions ("spin-dependent band 
shift" ) which guarantee the correct strong-coupling behaviour ||2ll . 

The present study focuses on the possibility and the stability of ferromagnetism 
in the PAM: Magnetic phase diagrams are constructed in terms of relevant model 
parameters such as hybridization strength V, /-level position Cf and total particle 
density n^^°^\ For this work we are mainly interested in the intermediate valence regime, 
i. e. that Cf is chosen within the energy region of the "free" Bloch band. The same holds 
for the chemical potential fi. As usual "/-level" and conduction band, respectively, are 
assumed as non-degenerate. The intraatomic Coulomb interaction U of the / electrons 
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leads to a splitting of Cf into two sublevels at Cf and Cf + f/. U is chosen so that the 
upper charge excitation remains unoccupied i. e. < n*^*"*-* < 3. The influence of V is 
multifold. First of all it provokes a finite width of the lower /-quasiparticle peak that 
turns out a basic prerequisite for ferromagnetism in the intermediate valence PAM. On 
the other hand, too strong s-f fluctuations due to V destabilize the ferromagnetic order. 
Additionally V provokes a hybridization gap in the energy spectrum, which, however, 
is getting closed by too strong quasiparticle damping. If the /-magnetization is almost 
saturated then the f-spectrum shows a hybridization gap, the J,-spectrum not. This has 
some influence on the induced s-polarization which can be parallel or antiparallel to the 
/-moment. The derived Curie temperatures are of reasonable orders of magnitude. All 
magnetic properties of PAM can be illustratively traced back and reasoned by inspection 
of the respective quasiparticle density of states. 
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